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ABSTRACT 

The Lick-index spectrophotometric system is investigated in its inherent statistical 
and operational properties to ease a more appropriate use for astrophysical studies. 
Non-Gaussian effects in the index standardization procedure suggest that a minimum 
S/N ratio has to be reached by spectral data, such as S/N> 5 px” 1 for a spectral 
resolution R ~ 2000. In addition, index (re-)definition in terms of narrow-band “color” 
should be preferred over the classical pseudo-equivalent width scheme. 

The overlapping wavelength range among different indices is also an issue, as it 
may lead the latter ones to correlate, beyond any strictly physical relationship. The 
nested configuration of the Fe5335, Fe5270 indices, and the so-called “Mg complex” 
(including Mgi, Mg 2 and Mgb) is analysed, in this context, by assessing the implied 
bias when joining entangled features into “global” diagnostic meta-indices, like the 
perused [MgFe] metallicity tracer. The perturbing effect of [OIII] 5 oo 7 and [NI] 5 ig 9 
forbidden gas emission on Fe5015 and Mgb absorption features is considered, and an 
updated correction scheme is proposed when using [OIIIJsocr as a proxy to appraise 
H/3 residual emission. When applied to present-day elliptical galaxy population, the 
revised H/? scale leads, on average, to 20-30% younger age estimates. 

Finally, the misleading role of the christening element in Lick-based chemical 
analyses is illustrated for the striking case of Fe4531. In fact, while Iron is nominally 
the main contributor to the observed feature in high-resolution spectra, we have shown 
that the Fe4531 index actually maximizes its responsiveness to Titanium abundance. 

Key words: techniques: spectroscopic - methods: observational - galaxies: elliptical 
and lenticular, cD - galaxies: photometry - galaxies: fundamental parameters 


1 INTRODUCTION 


Since its very beginning, the Lick-index system was inten¬ 
tionally conceived as a tool to c haracterize elliptical galaxies 
and other old stellar systems (Faber. Burstein. & Dressier 

- - -- "- Il98i 


197 


Faber et al 


1985i: 


Gorgas, Efstathiou. &; Aragon Salamancaf 


Davies et al. 

1199(3 ). 


Only 

a few relevant absorption lines were originally taken into 
account, at the sides of the strong Magnesium feature 
around 5175 A, which blends the molecular h ydride MgH 
and the atomic doublet Mgb ( |Mou 3 il978f) . Together 
with the 4000 A break, the Mg “valley” has since long 
been recognized as the most outstanding pattern i n the 
optical spectrum of galaxies (lOhma 3 11934 IWoodl 1 19631: 
ISpinrad fc Woodlll965 


* Also based on observations made at the Observatorio As¬ 
tronomico “G. Haro” of INAOE, Cananea (Mexico) 


The immediate interest for galaxy investigation di¬ 
rectly drove a parallel effort to extend the index database 
also to cool stars (e.g. iBurstein et alj fl984l ) for their evi¬ 
dent link with the study of old galaxies through popula¬ 
tion synthesis applications. At the same time, this also led 
to a more definite settlement of the whole index system, 
greatly extended by new additions from the Lick Image 
Dissector Scanner (IDS) data archive including all the rel¬ 
evant absorption features in the optical wavelength range 
comp rised between the 4000 A break and th e Balmer Ha 
line llGorgas et al.Hl9931 ; I Wort he v et al1ll994l ). The original 
emphasis to old stellar populations extend ed then to hot 
and warm stars (|Worthev fc Ottavianill 19971 ) to account for 
ongoi ng star formation scenarios, as in late-type galaxies 


miui iuimuixon ocoiiiM xuo, cvo ij.i ru,ro nj 

(e.g. Molla J= jIard^j_& = ^eauchampl!999 • Wal cher et al J[20()3 


iPerez, Sanchez-Blazauez. fc Zurital 20091 ). 

In more recent years, other spectroscopic methods have 
been complementing the original Lick approach, aimed at 
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tackling galaxy properties at a finer or coarser detail level. 
High-resolution spectroscopy (say at the typical resolving 
power, R > 2000 or so, of current extragalactic surveys) cer¬ 
tainly added valuable clues, on this line. On the opposite 
side, a panchromatic analysis relying on low-resolution fit¬ 
ting of galaxy SED (often assembled by straight flux con¬ 
version of integrated broad-band magnitudes) also became 
an increasingly popular way to deal with massive wealth of 
galaxy data, in order to characterize redshift distribution 
and other outstanding evolutionary features o n a statistical 
basis (see IWalcher et al.l l201ll : IConrovI 120131 . for excellent 
reviews on these subjects). 

For several reasons, however, none of these techniques 
is free from some limitations. From its side, high-resolution 
spectroscopy may be constrained by the intervening role of 
galaxy dynamics; at a typical velocity dispersion of a v ~ 
300 km s"\ the induced line broadening limits the effective 
resolving power at i?< 1000. On the other hand, global fit¬ 
ting of galaxy SED carries on all the inherent uncertainty in 
absolute flux calibration of the observations, especially when 
matching data from different observing sourc es and over dis¬ 
tant w avelength ranges. In addition, to quote IWalcher et al.l 
(1201 ill , we have to remind that, as “the uncertainties are 
dominated by the uncertainties in the SED modeling itself, 
thus one has to be very cautious about the interpretations 
when selecting samples where a specific type of model is 
preferred.” 

For all these reasons, narrow-band spectrophotometry 
still remains, in most cases, the elective (and often the only 
viable) tool to assess in finer detail the evolutionary status of 
galaxies and other unresolved stellar a ggregates both in lo¬ 
cal and high-redshift environments fseelBernardi et al . 20061: 
ICarson fc Nicholl2010l : iDobos et al.ll2012l . for recent contri- 
butions). Such a broad range of applications urges, however, 
a more careful in-depth analysis of the inherent properties of 
the Lick system in order to enable its appropriate use in the 
investigation of single stars and integrated stellar systems. 
This motivated the present work. 

When dealing with Lick indices, two intervening diffi¬ 
culties may in principle affect our conclusions. A first, and 
often underrated issue (discussed here in Sec. 2) deals with 
the index standardization , that is with the process to make 
index measurements, from different sources or observing cir¬ 
cumstances, to consistently compare eachother. Quite im¬ 
portantly, this process should not be confused with the cal¬ 
ibration procedure, as the first requires ab initio that index 
definition could lead anytime to a “fair” statistical realiza¬ 
tion of a normal statistical varia bl e acc ording to the Central 
Limit Theorem (see, e.g. ISac hsl l 19841, for a more rigorous 
theoretical settlement of the problem.) While a satisfactory 
attention is devoted in t he relevant literature to suitably 
calib r ate observing data (IWorthev et al.lll994l : iTraeer et al.l 
1 19981 : ICardiel et al.lll998h . this cannot assure, by itself, that 
a “standard” settlement of the derived index estimates is 
always achieved. 

The latter feature we have to remind is that indices 
may not be completely independent each one and, how we 
will show in Sec. 3, some degree of redundancy occurs in 
many cases. For instance, sometimes two indices may share 
the same pseudo-continuum window, or two feature windows 
are partly or fully superposed etc.. Even more sneaky, some 
redundancy may also appear when, under certain conditions, 


the same chemical element that constrains the main feature 
also affect the surrounding pseudo-continuum, thus masking 
or vanishing the index sensitivity to its abundance. When 
changing temperature or other fundamental parameters of 
stars, a “hidden” chemical element may abruptly appear at 
some point with its lines, and superpose to the other ab¬ 
sorption features affecting the corresponding indices. This 
is, for instance, the classical case of the onset of TiO molec¬ 
ular bands in cool sta rs perturbi ng , among others, t he 


Mg 2 a nd Fe5270 indices (BuzzonL_Gariboldi, & Mantcgazza 
1992|: Fhantalo fe Chiosil 20041 : iBuzzoni. Bertone. fe Chavez 


2009all . 


A further range of problems deal with more physical 
environment conditions of stellar populations. Most impor¬ 
tantly, a natural effect of fresh star formation on the in¬ 
tegrated spectrum of galaxies is line emission driven by the 
residual interstellar gas. Along the Lick-index range, this re¬ 
sults in enhanced Hydrogen Balmer lines and, depending on 
the overall thermodynamical conditions, also in [ 0111] for- 
bidden emission at 4363, 4959 and 5007 A (see iKenniceittl 
119921. for excellent explanatory examples). On the other 
hand, Oxygen and H(3 emission is triggered even at the van- 
ishi ng star formatio n rates observed in m o st elliptical galax- 
ies I Caldwell 19841 : Phillips et al.l 19861 : Bettoni fc Busoni 
1 1987 ; Volkovl Il990l : ISarzi et ahl 20ld : Crocker et al.1 201lh 
due, at the very l east, to the contribution of the planetary 
nebulae fsee lKehrig et al.ll2012l . for an example), and this 
could easily bias our conclusions if we use H/3 (absorption) 
strength to constrain galaxy age. 

In this paper we will carry on our analysis from two 
independent points of view, starting from the expected be¬ 
haviour according to index properties and checking theoret¬ 
ical predictions by means of a tuned set of spectroscopic 
observations for a sample of twenty elliptical galaxies. Our 
relevant conclusions, and a reasoned summary of the the 
manifold variables at work is proposed in Sec. 4 of the pa¬ 
per. 


2 INDEX TECHNICALITIES 


Up to the 70’s, low-resolution sampling of galaxy and stel¬ 
lar energy distribution was mainly pursued by narrow- 
band photometry. P opular mul ticolor systems included, 
fo r instance, those o f Istrorng ren] (|l956l h the DDO system 
oflMcClure fc van den Berghl ( 1968l l. the multicolor set of 
S pinrad fc TavloiT i 19691 )' or the 10-color system of iFaberl 
lll973ll . In this framework, the Lick system stands out for 
its innovative approach to the problem, tackled now from a 
purely spectroscopic point of view. 

Indeed, such a new “spectroscopic mood” is inherent in 
the Lick-index definition. The procedure applies in fact to 
low-resolution (R ~ 500) spectra and aims at deriving the 
strength of a number of absorption features through a mea¬ 
sure of their pseudo-equivalent width. Indices are therefore 
measured in A, maintaining the ma gnitu de scale only for 
those few cases of molecular bands (IWorthev et alj 199411. 
The original set of 21 indices was defined by Worthev et al.1 
(Il994ll and slightly a dapted to mo r e suit ably fit with ex¬ 
tragalactic targets by iTrager et al.1 (119981 4 . A complement 
to include the full Bal mer li ne series i n the optical wave¬ 
length range led IWorthev fc Ottavianil (Il997l f to add four 
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Figure 1. A global vie w of the Lick-index wavelen gth distribution. The original set of 21 indices dWorthev et al.lll994f) has been 
comple mented by the fou r lWorthev Sz Ottavia ni (1997) new Balmer indices (namely, H5 and H 7 , in their “wide” and “narrow” version) 
and the lGonzale3 (1 19931) [OIII] index, as labelled in the plot. The i ndex map is superposed to a synthetic stellar spectrum with 
[T e ff, log g, Z] = [4750 K, 2.00, Zq] from the Bluered theoretical library llBertone et a.l.[12004! . 12008 ). degraded to match the Lick spectral 
resolution of 8.5 A FWHM. 


new indices, that match H'y and HS lines with “narrow” and 
“wide” windows. With these additions, the full Lick system 
consists of 25 indices. Apart, it is worth including in our 
analys is also an extem porary addition to the standard set 
due to lGonzaled (1 19931 ). in order to account for [OIII] emis¬ 
sion in the spectrum of star-forming galaxies. Note that this 
is the only index aimed at measuring an emission feature 
(see Fig. [[]). 

According to IWorthev’sl (Il994l ) original precepts, the 
key quantity of any Lick index is the ratio 



1Z is therefore an average, along the w-wide feature win¬ 
dow, of the apparent flux density, /(A), normalized to the 
(pseudo)-continuum level, / c (A), linearly interpolated at the 
same wavelength. Two adjacent bands, on both side of the 
feature, provide the baseline for the interpolated process. 
With these prescriptions, the index then directly follows ei- 


ther as 


IT' 

1 

l—1 

II 

°5l 

(2) 

or 


/(mag) = —2.5 log 1Z 

(3) 


depending whether we want to express it in A, as pseudo¬ 
equivalent width, or in magnitude scale. Clearly, both defi¬ 
nitions are fully equivalent and easily reversible as 


/(A) = W [l - io-°- 41(mag) ] 

(4) 

or 


/(mag) = -2.5 log [l - -i^l, 

L w J 

(5) 


and the use of A or magnitudes is just a matter of 


“flavor”, as we have been men t ioning above (see, for 
instance, Brodje_J^JIuclira 1990l . IColless et al.l 1 19991 . or 
ISanchez-Blazauez et al. 2006a for alternative non-standard 
notations). As a rule, both eq. @ and 0 are positive for 
absorption features and negative for emission ones. 


2.1 Standardization constraints 


In spite of its established popularity and extensive use 
in a broad range of stellar and extra-galactic studies, 
a well recognized drawback of the Lick system is its 
ambigous standardization procedure. A meticul ous an al¬ 
ysis i n this sense has been attempted by iTrager et alj 
(| 1998t ) to account for the macroscopic limits of the IDS 
database due to the somewhat fickle technical performance 
of the instrument. Some inconstancy in the spectrum 
wavelength scale is, for instance, a recognized issue, and 
even th e spectr al resolution (about 8.5 A FWHM, see 
iWorthev fe Ottavianil 1 19971 ) cannot be firmly established 
due to slightly un predictable behaviour of the instrument 
output with time (IWorthev et al.l Il994l ). In addition, and 
quite importantly, one has to notice that t he wh ole IDS 
data set is in fact not flux-calibrated (IWorthev et alJ 
1199dl) . thus taking over the instrument (slightly variable) 
response curve. For all this, any comparison between 
external data sources in the literature must forcedly rely, 
as a minimum, on th e observation of a common grid of 
stand ard s tars (e.g. Worthe^et^L _ 1994 1 _ Gorgas et al.1 
1993|j _ Buzzoni, Gariboldi. fe Mantegazza 1992 ; 


Buzzoni. Mantegazza. fe Gariboldil 1 1994 [b 

2001, among others). 


li et al. 


A few attempts have been made, in the recent literature, 
to refine the index estimates from the operational point of 
view to secure on mor e solid bases the observed output. The 
work of lRogers et al.l (l2010i) is certainly the most relevant 
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one, in this sense, as it especially focus on the problem of 
a fair settlment of the pseudo-continuum level (through the 
so-called “boosted median continuum” method) to derive a 
more confident feature strength, also overcoming those rec¬ 
ognized cases of unphysical output (for example with nomi¬ 
nally “emission” indices tracing, in fact, absoption features 
etc.). 

As a matter of fact, any more or less elaborated empir¬ 
ical matching procedure, however, cannot overcome a more 
subtle but inherent difficulty in the standardization process 
of the Lick system. In fact, the way one leads to compute 7Z 
according to eq. © is a direct “heritage” of high-resolution 
spectroscopy, where the measure is always carried out on 
preliminarily continuum-rectified spectra. However, this is 
clearly not the case for the (unfluxed) Lick spectra, where 
possibly steep continuum gradients are in principle an is¬ 
sue. According to statistics theory, the consequence is that 
ratio 1Z may d isplay strongly non-Gaussian properties (see 
lMarsiglialll964l . and especially Fig. 1 therein for an eloquent 
example) making the analysis of its distribution a hard (and 
in some cases hopeless) task. In other words, this implies 
that, in general, any effort to reproduce the system may not 
obey the Central Li mit Theorem so that its convergence is 
not assured a priori (jMarsiglialll965l l. 

In this framework, it is therefore of central interest to 
explore the real contraints that make 7 Z behave normally. 
For this, let us briefly consider the more general problem of 
the stochastic behaviour of a statistical variable r being the 
ratio of two normal variables x and y , nam e ly r = y/x. It 
has been firmly demonstrated (IGearvI Il930l : iHinklevI 19691 : 
IHavva. Armstrong, fe Gressislll975ll that r tends to behave 
like y if the range of variation of x (that is dx/x), approaches 
zero and the x va riable is “unlikely to assume negative val¬ 
ues” llGearvil93(tl . If this is the case, then x tends to behave 
as a (posi tive) “constant” and xr wou ld resem ble the y dis¬ 
tribution dHavva, Armstrong, fc Gressislll975l h Under these 
assumptions, the r' = ylx ratio will display a normal distri¬ 
bution and its average r' tends to be an unbiased proxy of r. 
In our terms, therefore, the crucial constraints that ensure 
both index Gaussianity and its fair reproducibility directly 
deals with two issues: from one hand we need a “small” range 
of variation for df c / f c ; on the other hand, we need to assess 
the condition for which r' —I r. 

To further proceed with our analysis, let us preliminar¬ 
ily assume the relative variation of f c to be “small”, indeed, 
so that we are allowed to set up a linear expansion for f c 
around A 0 , the central wavelength within the feature window 
of width w: 


1 

To 



1 



dfc 

d\ 


(A — A 0 )] 


1 

% 



a{X _ K) ]. ( 6 ) 


In previous equation, a = df c /d A and, by definition, (aw) = 
df c is the maximum excursion of f c ( A) within the feature 
window. 

By implementing eq. © into eq. © we have 


n= lrUWdx 

fc L w 


4= f /(A)(A-A„)dAl. 

wfc J w ^ 


(7) 


If we now multiply and divide by w f f(X)dX the last term 
of the equation, with little arithmetic we can re-arrange the 
equation and write 


^ aw 

& 

II 

i-M, 

/ w \ max 

i) 1 

^ fc 

W J L 

w 

\fc) J 


being 


(8) 


tt'=(=), (9) 

and (SX)/w the relative displacement, within the feature 
window, of the line centroid with respect to A 0 , in conse¬ 
quence of symmetry departure due to the spectral slope. 
According to eq. ®, therefore, 1Z' tends to become a fair 
proxy of 1Z if 




( 10 ) 


The maximim relative excursion (d/c//c) max is, evi¬ 
dently, the key figure in this process as it directly constrains 
both r.li. factors in eq. m- It may be wise to assess the 
allowed range of this parameter in terms of the stochastic 
fluctuation range of f c within the same window, the latter 
being directly related to the signal-to-noise ratio ( S/N ) of 
our data, so that we can set 



with k a tuning factor. We will show, in the next section, that 
a minimum signal-to-n oise ratio o f S/N ~ 20 is requir ed for 
7 Z! to behave normally dHavva, Armstrong, fc Gressidll975h . 
With this figure, if we choose to tentatively accept a pseu¬ 
docontinuum excursion within a ±3<r range of the random 
noise, then fc< 6 and (d/ c // c ) max < 0.3, so that one could 
actually verify numerically that 


1 - 


TZ_ 

7^' 


< 0 . 01 . 


( 12 ) 


With the previous assumptions, this means that, within a 
1% accuracy, the ratio 1Z' eventually becomes an unbiased 
estimator of 1Z, thus assuring 1Z to beh ave n or mally , too. 
Our first conclusion is therefore that I Gearv ’ j (1 193G ) pre¬ 
scriptions require pseudo-continumm not to exceed a ±30% 
relative variation within the feature window for any Lick in¬ 
dex to obey the Central Limit Theorem, and display a Gaus¬ 
sian distribution. 


2.2 Random-error estimate 


A range of more or less entangled attempts have been pro¬ 
posed across the literature to derive a general and self- 
consistent estimate of the sta tistical uncertainty to be re- 
lated t o narrow-ba nd in dices dRichl 1 19881: iBrodie fc Huchnil 


jI990|^_ Carollo 1 J3anzj^r 1 _& L ^^on|_ 19931 : Gonzalez 1993; 
Cardiel et al.ll 19981 : iTrager et al.lll998fl . Yet, this effort does 
not seem to have led to any firm and straightforward path¬ 
way to assess the index random-error estimate. Quite differ¬ 
ent formal approaches can be found among authors depend¬ 
ing whether the intervening noise sources (including photon 
statistics, flat-fielding procedure, sky subtraction etc.) are 
considered individually or to a somewhat aggregated level 
in the analysis. In addition, as we have been discussing in 
previous section, with the Lick system one has to deal with 
the duality of index definition, sometimes seen in terms of 
pseudo-colors and sometimes in terms of pseudo-equivalent 
widths. 
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In its general form, by simple differentiation of eq. m 
and m, the expected internal uncertainty of a Lick index 
expressed in magnitude is 



On the contrary, if we chose to express the index in pseudo¬ 
equivalent width, then 


< 14 > 

In both cases, the crucial quantity that we have to as¬ 
sess is the o(JZ)/7Z ratio. A usef ul contribution t o settl e 
the problem has come from IVollmann fe Eversberd l|2006l ) , 
dealing with the error estimate in equivalent-width measure¬ 
ments from high-resolution spectroscopy. Under many as¬ 
pects, their results directly apply to our framework provid¬ 
ing the caveats of previous section to be properly accounted 
for in our analys is and to observe at “sufficie nt S/N within 
the continuum” llVollmann fe Eversberdfeood f. If the 7Z 1 ra¬ 
tio could be an effective proxy of 7Z, then the variance of 1Z 
directly derives from error propagation as 

<»> 

which implies 


- (f/i) + ( 



) 


f2 2 n2 ^-2 

_ ZL i 2 _ 

f 2 f2 ' f2 f 2 ' 

J c J J c J c 


(16) 


In case a similar S/N ratio along the relevant range of our 
spectrum can be assumed, so that (cr(/)//) ~ (ct(/c)// c ), a 
final compact form can be achieved for previous equations 
such as 


u(JZ) ^ &{TZ) \/2 

~lZ~ ~ ~W = (S/N) idx ■ 


Quite importantly, note that the value of S/N in 
eq. HD refers to the integrated figure along the feature 
and continuum windows. Providing to observe a spectrum 
with an original ( S/N) 0 b s ratio per pixel and a dispersion 
8 A px _1 , the resulting ( S/N)i dx by sampling flux along a 
window w A wide is, of course, 


{S/N)id x = ( S/N ) obs ./H (18) 

V & 

In their work, dHavva. Armstrong, fe Gressislll975ll set 
the limit of general validity of this formal approach through 
a detailed Monte Carlo analysis. Converting their results to 
our specific framework (see, in particular, their Fig. 1), we 
conclude that eq. El) is a robust estimator of the relative 
standard deviation of 7Z providing to work with ( S/N)i dx > 
20 . 

On the basis of our discussion we can eventually re¬ 
arrange previous eq. m and El in their final form, re¬ 
spectively: 

&(I) mag 
■ 

m = [?fL w 

a[ ’a Vw (S'/A/Us 
For “shallow” lines (that is those with [7 (A)/ui] -C 1), we 


^ ■ 


( 20 ) 


1.09 


(S/N) 0 


(19) 


Table 1. Harmonic average^) (W), feature (w) and continuum 
(W) window width for the extended set of Lick indices 


Index W w W Index W w W 

- [A] - - [A] - 


H(5yi/ 

50.1 

38.8 

70.6 

H5jv 

30.5 

21.3 

53.8 

CNi 

48.2 

35.0 

77.5 

cn 2 

42.0 

35.0 

52.5 

Ca4227 

15.0 

12.5 

18.8 

G4300 

33.7 

35.0 

32.5 

H 7w 

58.6 

43.8 

88.8 

Hqjv 

31.9 

21.0 

66.3 

Fe4383 

32.5 

51.3 

23.8 

Ca4455 

23.1 

22.5 

23.7 

Fe 45 3i 

35.1 

45.0 

28.8 

Fe 46 68 

47.3 

86.3 

32.5 

H/3 

31.6 

28.8 

35.0 


[OIII] 

32.7 

20.0 

90.0 

Fesois 

54.6 

76.3 

42.5 

Mgi 

86.1 

65.0 

127.5 

Mg 2 

63.8 

42.5 

127.5 

Mgi, 

33.1 

32.5 

33.7 

Fe5270 

43.4 

40.0 

47.5 

Fe5335 

27.8 

40.0 

21.2 

Fe5406 

24.0 

27.5 

21.3 

Fe5709 

29.1 

23.8 

37.5 

Fe5782 

21.7 

20.0 

23.7 

NaD 

36.3 

32.5 

41.0 

TiOr 

72.3 

57.5 

97.5 

Ti0 2 

97.0 

82.5 

117.5 


According to eq. GD- 


can simply neglect the last factor in previous equation, thus 
leading to an even more compact (though slightly overesti¬ 
mated) form for the standard deviation. 

In both eq. El and (pO)). W is the harmonic average 
of the feature (w) and continuum (W) window widths: 


W “ \w + w)' 


( 21 ) 


For reader’s better convenience, the value of W, together 
with the feature and the full (i.e. red+blue) continuum win¬ 
dows for the extended set of Lick indices is summarized in 
Table U 

Taking eq. El and (1201) as a reference, in Fig. [2] we 
compare the observed standard deviatio n of pr i mary Lick 
stellar calibrators, according to IWorthev et al. d 19941. see, 
in particular, Table 1 therein) and Worthev fe Ottavianil 
(|l997l ) with our expected error distribution for a spectrum 


of (S/A0obs = 30. 

Although no firm appraisal can be done of original IDS 
data in terms o f S/N ratio, a problem extensively discussed 
bv IWorthev et all (Il994lh nonetheless there is point to be¬ 
lieve that the S/N assumed in the present derivation is sim¬ 
ilar to the floor imposed on the IDS spectra by small-scale 
flat-fielding problems. In this framework one could recog¬ 
nize our statistical analysis as an instructive excercise to 
put in a more standard context the inherent uncertainty of 
Worthey’s et al reference dataset, at least as far as indices 
“redder” than ~ 4350 A are concerned. 

Note, however, a dramatic worsening of the IDS perfor¬ 
mance in the blue, where index accuracy quickly degrades. 
This trend can now easily be understood on the basis of 
our previous statistical arguments. At shorter wavelength, 
in fact, the combined effects of a worse IDS response and an 
intrinsic drop of SED for stars of intermediate and late spec¬ 
tral type all conspire in the sense of exacerbating any spec¬ 
tral slope, thus departing from the standard Gaussian sce¬ 
nario (in force of eq. El, and adding further extra-volatility 
to Lick-index calculation. 
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Figure 2. The expected index standard deviation cr(I) 0 ur , ac¬ 
cording to eq. Eoi . for an assumed spectral S/N ~ 30 is com¬ 
pared with the empiri cal estima t e from the Lick primary calibra¬ 
tors, a s from Tabic 1 of Worthcy (1994) and Wor thev &: Ottavianil 
d 1997ll - The uncerta inty for [OIII] 5 oo 7 has been estimated from 
the iGonzale j (Il993l l data, as from Table 4.1 therein, by scaling 
down his figures to S/N = 30. Notice that Lick IDS spectra lack 
any stable intensity response, a feature that pre vents any firm as¬ 
sessment of the S/N ratio o f the original data iFaber fe Jackson! 
1976: Wor thev et all 1 1 9941 ) . A remarkably good agreement is 
found, however, with our theoretical predictions, as far as “red” 
indices (about A >4350 A) are concerned. A notable worsening 
in the Lick calibration accuracy is evident at shorter wavelength, 
likely due to poorer detector performance and statistical draw¬ 
backs inherent to index definition, as discussed in the text. 


3 INDEX ENTANGLEMENT AND 

REDUNDANCY 

A closer look to Fig. |T] data shows that about 25% of the 
covered wavelength range actually contributes to the defi¬ 
nition of two or more indices leading, on average, to a 1.5 
oversampling fraction. Such a redundancy evidently calls for 
a more detailed analysis of the possible bias effects that in 
some cases make indices inter-dependent in their stochastic 
fluctuation, that is beyond any strictly physical relationship. 

This is certainly the case of some intervening emission, 
that superposes to absorption features ( H/3 is a typical ex¬ 
ample in this sense), but also partly or fully overlapping 
pseudo-continuum windows, in common with different in¬ 
dices, may induce some degree of correlation in the corre¬ 
sponding output (as for the blue sidebands of Ca4455 or 
Fe5335, fully comprised within the Fe4383 and Fe5270 red 
sidebands, respectively). This may also happen to feature 
windows, as shown for instance by Mgb, fully comprised 


within the Mg 2 central band. Finally, a further possibility is 
that the entire feature contributes to the pseudo-continuum 
window of another index: the case of Fe5335, is explanatory 
in this regard being fully nested into the red sideband of 
Mg 2 and Mgi. 

In this section we want therefore to assess some illus¬ 
trative cases giving, when possible, a guideline for a more 
general application of our results. The theoretical predic¬ 
tions that come from the reference discussion of Sec. 2, will 
be accompanied by an empirical check, that could allow us a 
clean assessment of the different situations sketched above. 


3.1 Radial spectroscopy of elliptical galaxies as a 
robust empirical check 

For our aims it is convenient to rely on a set of spectro¬ 
scopic data for a sample of 20 early-type (mainly elliticals) 
galaxies observed with the 2.12m telescope of the “G. Haro” 
Observatory of Cananea (Mexico) in a series of runs along 
years 1996/97. These observations are part of a more gen¬ 
eral exhaustive study on spectroscopic gradients across the 
surface of gala xies and its physical implications for galax y 
evolution (ICarrasco et aljfl995l : iBuzzoni et al.|[2009bl . 20151 ). 
It is clearly beyond the scope of this work to dig further into 
the original scientific issue and the detailed data reduction 
analysis (the reader could be addressed to the cited refer¬ 
ences for a full review). We just want to recall here the basic 
information to characterize the sample in view of its possi¬ 
ble relevance as an independent tool to probe our theoretical 
output. 

The spectra have been collected with a long-slit Boiler 
& Chivens Cassegrain spectrograph equipped with a 300 
gr/mm grating, that provided a dispersion of 67 A/mm 
along a 4200-6000 A wavelength interval. Spectral resolu¬ 
tion was set to 5 A (FWHM) throughout (R = A/A A ~ 
1000) Galaxy targets were typically exposed along one hour 
(3 x 1200 ,/ frames) centering along galaxy major axis. The 
resulting spectra have been wavelength- and flux-calibrated 
according to standard procedures, and have then been re¬ 
binned in the spatial domain such as to sample galaxy ra¬ 
dius at 3-5 values of distance symmetrically placed with 
resp ect to the center, ro ughly up to one effective radius 
(see lBuzzoni et ail l2009ll . for an example). The rebinning 
procedure allowed us to greatly improve the S/N ratio es¬ 
pecially in the outermost regions such as to work there 
with S/N >25 px -1 , a figure that quickly improved to 
S/N >80 px -1 toward the center. 

The relevant information of the galaxy sample 
is summarized in Table [5] From the RC3 catalog 
(Ide Vaucouleurs et al. IH) we report morphological clas¬ 
sification, absolute B-band magnitude (Ms), and apparent 
(B — V ) color within one effective radius (r e ). In addition, 
in the last two columns we added, respecti vely, the radial 
fractio n sampled with our spectra (r/r e ) (see lBuzzoni et al.l 
l 2009bl . for further details), and a flag for galaxies with re¬ 
ported emission. In this regard, o ne asterisk mark a con - 
firmed [ 011)3727 emission, as in the lBettoni fe Busoni (Il987l) 
catalog, while with a double asterisk we marked targets with 
evident [OIII] 5 oo 7 and/or H/3 ce ntral emission in our spectra . 

For the galaxies of Tablel2l lBuzzoni et al] (l2009bl . f2015h 
have derived the complete set of 17 standard Lick indices 
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Table 2. The galaxy database 


NGC 

Type 

M s 

(B-V) 

r max Ae 

emission? 

1587 

Elp 

-21.19 

1.04 

0.2 

* * 

1588 

Ep 

-19.91 

1.00 

0.5 


2685 

SO-a 

-19.09 

0.94 

0.8 

* 

2764 

SO 

-19.55 

0.71 

1.2 

* * 

3245 

SO 

-20.09 

0.89 

1.0 


3489 

SO-a 

-19.23 

0.85 

0.8 

* * 

3607 

E-S0 

-20.00 

0.93 

0.4 

* 

4111 

SO-a 

-19.16 

0.88 

3.0 

* 

4125 

E6p 

-21.27 

0.94 

0.4 

* 

4278 

El 

-19.35 

0.96 

1.0 

* * 

4374 

El 

-21.05 

1.00 

0.6 

* 

4382 

SO-a 

-20.50 

0.89 

0.3 


4472 

E2 

-21.71 

0.98 

0.8 


4649 

E2 

-21.47 

1.00 

0.5 


4742 

E4 

-19.35 

0.79 

1.3 

* 

5576 

E3 

-20.14 

0.90 

0.4 


5812 

E0 

-20.44 

1.04 

0.4 


5866 

SO-a 

-20.00 

0.93 

0.7 

* * 

5982 

E3 

-21.38 

0.94 

0.7 

* 

6166 

E2p 

-22.96 

1.03 

0.4 

* 


(including the [OIII] 5 oo 7 feature), from G4300 to Fe5709. 
These data provide us with a useful reference to perform, 
for each index, a series of statistically equivalent realizations. 
This can be carried out in force of the supposed central 
symmetry of the galaxy spectral properties. Even in case 
of spectroscopic radial gradients, in fact, a similar stellar 
population has to be hosted at the same distance on opposite 
sides of the galaxy leading to a nominally identical set of 
absorption indices. Clearly, any systematic deviation from 
this pure folding effect may be worth of attention in our 
analysis. 

Note that, as we deal with index differences, our method 
is quite insensitive to standardization problems (zero points 
settlement, calibration drifts etc.), although we obviously 
magnify the internal random errors of our analysis by a fac¬ 
tor of y/2. About 80 couples of index differences have been 
computed across the full galaxy sample, for a total of roughly 
1400 measures for the full set of 17 indices. 

Just as an explanatory example, to illustrate our ap¬ 
proach, we show in Fig. [3] the A-A index distribution for 
two Fel indices, namely Fe5270 versus Fe4383. Each point 
in the plot comes from the following procedure. For two 
slices at the same radial distance (on opposite side) of a 
given galaxy, we computed 8x as the “left - right” differ¬ 
ence of the Fe4383 index measurements. The same is done 
for Fe5270, deriving 5y. We then set Ax = |5a:| and obtain 
the sign of 5x as sign^. = 5(x) / A(x), so that Ay = signal/. 
Though, from a physical point of view, both indices should 
correlate with Iron abundance, no mutual influence can ex¬ 
ist for their relative variation. As a consequence, their “left- 
right” differences must be statistically independent over the 
entire galaxy sample. This is actually confirmed by the least- 
square fit to the data in the figure (68 points in total), which 
provides a = 0.03±o. 05 0 


o 

(M 

iO 

< 1 ! 

fc, 

< 


0 2 4 6 8 

AFe4383 

Figure 3. T he envisaged A-A test amo ng galaxy spectroscopic 
indices in the (jBuzzoni et all2009bl . l2015h galaxy sample (see text 
for details) is explained here for the illustrative case of Fe4383 (in¬ 
dex “x”) and Fe5270 (index “y”). If we assume the source of “left 
- right” index difference to be only of stochastic nature, then no 
correlation has to be expected between A Fe4383 and A Fe5270. 
A least-square linear fit to the data, in the form of Ay = a Ax, 
actually provides a = 0.03±o.O5 (solid line with the yellow fan 
marking the drier standard deviation of the slope estimator), fully 
consistent with a flat trend (dashed line), as expected for this case. 
Just as a guideline, notice that the dot size in the plot is smaller 
for index measurements at farther radial distances from galaxy 
center. 


“i-1-1-1-1-1-1-1-1-1-1-1-1-1-r 


teo: AFe52= 0.00 AFe43 “ 



• • 


3.2 Menage a trois: H/d, [OIII] 5 oo 7 and Fe5015 

The wavelength region of galaxy spectra around 5000 A may 
be heavily prone to the effects of gas emission. In particular, 
the H/3 and Fe5015 strengths could be affected, the first 
resulting from the parallel action of stellar absorption and 
gas emission, the latter because of the possible [OIII] 5007A 
appearence (see Fig. [4j) . This issue is of central relevance as, 
within the Lick system, the HP index is the most suitable 


stellar aggregates (e.g. 

Buzzoni. Mantegazza. & Gariboldi 

1994 Bressan, Chiosi, & Tantalo 199q; Vazdekis et al. 1996|; 

Tantalo. Chiosi. & Bressan 19981; Kobavashi & Arimoto 

1999 

; Maraston & Thomad 2000; Trager et al. 

200d 

: iThomas et al. 

20051; Fusi Pecci et al. 20051; 

Sanchez-Blazauez et al. 

2006bi: Peletier et alJ 200?l. and 


many others) due to its close response to the temperature 
locat i on of the Main Sequence Turn Off stars (|Gorgas et al.l 
1 19931: IB uzzon ] |l 9951 : [rUizzoni. Bertone. &; Chavedl2009al f~ 


A proper correction of this effect is not quite a straight¬ 
forward task, especially when modest amounts of emission 
do not appear as “surging” features in the spectra and 
rather hide making absorption line-depth artificially shal¬ 
lower. If uncorrected, this effect would bias age estimates 
toward older values. A fir st attempt to ac count for H/3 e 
residual emission is due to I Gonzalezl dl993l l , who assumed 
the [OIII] 5007 index to be a confident proxy for Hj3 cor- 


the catastrophic outliers (residual spikes from cosmic rays in the 
1 As a general rule, throughout in our A-A fits, we applied both to spectra, cold/warm pixels affecting the index strength, lacking 

x and y variables an rms clipping procedure at a 3i r level to reject index estimate on either “left” or “right” side etc.). 
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Wavelength [A] 

Figure 4. The nested configuration of the [OIII ]5007 and 
Fe5015 indices. For each index the two blue and red pseudo¬ 
continuum sideband are reported, flanking the main feature win¬ 
dow. Two representative spectra of active (NGC 2764) and quies- 
cent (NGC 4649) ellip tical galaxies, from our observed database 
(IBuzzoni et alj l2009bh , are superposed to the index sketch, for 
the sake of comparison. In addition to the [OIII] emission, note 
for NGC 2764 also a prominent Hemission at the nominal 
AA 4861 A. 


rection. A simple scheme was devised in this sense, where 
the observed Balmer i ndex sho uld be offset by AH/3 = 
Hf3 e = — 0.7 [OIII] 5007 llGonzalejll993h . This approach re¬ 
lies on the fact that [0111] emission would be much more 
clearly recognizable than possibly hidden H[5 emission, and 
[OIII] and Balmer emission do somewhat correlate as they 
ostensibly originate in the same regions of the galaxies. 
One has to note, howe ver, that later attempts to settle 
the problem empirically llCarrasco et al.l Il995l : il Yager et al.l 
120001 ; Kuntschner et al .1120061 ) reached puzzling conclusions, 
at odds with Gonzalez’ predictions. 

To a closer analysis, things may actually be much 
more entangled as i) apparent [OIII] emission is partly 
at odds with Fe5015 absorption, as both indices partly 
overlap (see Fig. [4}: ii) fresh star formation could likely 
trigger both [OIII] _and Bal mer emissio n b ut o n differ- 
ent timescales llMoustakas. Kennicutt. fc Tremontil 120061 : 
iBuzzoni et al .1 12 01 5l ). and de pending on Oxygen abundance 
and i onization parameters (lOsterbrockl 1 19741 : IPagel et al.l 
Il979l : iKewlev fc Dooital 12002 1. For all these reasons, the 
combined behaviour of the three indices, H/3, [OIII] 5 oo 7 and 
Fe5015, needs to be assessed self-consistently. 

It is useful to start our analysis first with the [OIII] 5007 
versus Fe5015 interaction. Figure [4] clearly shows that Oxy¬ 
gen emission may be systematically underestimated as the 
[OIII] 5007 continuum sidebands cannot account for Fe5015 
absorption. This is especially true if o ne considers that 
Fe5015 strengthens in old stellar system s (|Tager_et_alJ|2000|; 
iBeuing et, allhooil lOgando et al.lf2008l : IBuzzoni et aljhoifih 
(see, for instance, the illustrative case of galaxy NGC 4649, 
also displayed in Fig. 0. Let us then compute the observed 
flux density within the [OIII] and Fe5015 feature windows 
(of width wo and W50, respectively). 

By recalling eq. m and m, and omitting the “aver- 



Figure 5. The A-A correlation test, as for Fig. [3] for the 
[OIII ]5007 and Fe5015 indices. The expected statistical prediction, 
according to eq. J28D (dashed line), is well matched by the obser¬ 
vations, the latter ones providing a fully consistent a = 1.15±o.2i 
fitting line slope (solid line within a Tier yellow fan), as reported 
top right in the plot. 


age” overlined notation for flux density symbols, as well as 
the suffix to the Oxygen index, for better legibility of the 
formulae, we have 


fo_ _ A _ /so A .Fe5015 \ 

fc V wo / ’ fc \ W50 '' 


( 22 ) 


Note that fo is in fact a sum of the genuine contribution 
of O emission plus the contribution of the intrinsic Fe5015 
feature. Similarly, the intrinsic Fe5015 flux density (i.e. after 
fully removing O emission) is 


/50 = A Fe5015' 

fc v W50 


(23) 


where Fe5015’ stands for the corrected Fe5015 index 
strength. By definition, it must be 


/5o(w50 — wo) = /50 W 50 - fo wo 
so that, after matching eq. <122)1 . 


/50 

fc 


(wso — Wo) 


w O ) 


(24) 


(25) 


When combining eq. (1231) and os this lead to 

Te5015' = W 50 = ——- (Fe5015—[OI 11]) .(26) 

V fc / w 50 — wo k V 


By replacing the corresponding width of each feature win¬ 
dow (namely wo = 20 A and wso = 76.3 A, see TableQ]), the 
intrinsic Fe5015 index can eventually be written explicitely 
in terms of the observed quantities as 


Fe5015' = 1.36 (Fe5015 - [O///]) . (27) 

As far as the A-A index distribution is concerned, if the 
same (intrinsic) value of Fe5015’ has to be expected in the 
galaxy “left” and “right” side, at fixed radial distances, then 
any difference in the observed Fe5015 index may be induced 
by a change in the [OIII] emission. By simply differentiating 
eq. m, we obtain a straight positive correlation such as 
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Fe5015 [A] 

Figure 6. The ob served (raw) [OIIII 50 07 versus Fe5 015 index 
distrib ution for the iQgando et alJ ( l2008h (s olid dots), lGonzale3 
(1 19931) (triangles) and (iBuzzoni et al. II 200913 , 120151) (rombs) sam¬ 
ples of early-type galaxies. The forbidden zone in the index do¬ 
main is shaded, according to eq. (3D, while a more general pa¬ 
rameterization versus intrinsic Fe5015 and [OIII ]5007 strength 
(primed indices in the plot) is displayed, in force of eq. (3D and 
ED. An indicative error box for the observations is reported top 
right. 

AEe5015 = A [Oil I]. (28) 

Again, when comparing in Fig. [5] with our data, a least- 
square solution on 68 selected points in total (after a 3cr clip¬ 
ping) provides a fully consistent slope value a = 1.15±o.2i- 
Quite importantly, turning back for a moment to 
eq. (1271) . one has to notice that, when [OIII] emission fades 
and Fe5015 tends to Fe5015’, then the [OIII] index does not 
vanish. In fact, if we set .Fe5015 = Fe5015' in eq. Hi , one 
obtains 

[OIII] min = —Fe 5015' = 0.26 Ee5015' (29) 

who 

More generally, if some emission occurs, then [OIII]' = 
[OIII] - [0///] min , and 

[OIII]' = [OIII] - 0.26 Fe 5015', (30) 

or 

[OIII]' = 1.36 [OIII] - 0.36 Ee5015 (31) 

putting it in terms of the observed quantities. Equation Hi 
gives the locus in the Fe5015 versus [OIII] plane for any 
intrinsic emission strength [OIII]'. 

By combining eq. Hi and pj . an obvious self- 
consistency constraint implies that 

Ee5015' = Ee5015 - [OIII]'. (32) 

A graphical summary of our theoretical scheme is 
sketched in Fig. [D] where we compare the [OIII] versus 
Fe5015 distribution for the galaxies of Table [2] with the ob- 
serve d distribution o f an e xtended sampl e of 509 ell i pticals 
from IQgando et al.l (l2008l l and with the iGonzalexl (Il993l l 
original bulk of 41 ellipticals. The trend seems quite en¬ 
couraging, and actually calls for an important property of 



1 0-1-2 -3 


[OIHl'aoo, [A] 

Figure 7. The [QUIT 5 007 distribution of the intrinsic (“primed”) 
index among the IQgando et alJ d20Q8h early-type galaxy sample. 
A total of 498 galaxies, with available [OIII] 5007 and Fe5015 ob¬ 
servations are displayed. The [OIII] observations have been cor¬ 
rected according to eq. ED- On average, ellipticals display a weak 
residual gas emission with {[ OIII ]' 5007 ) = — 0.70±o.54 A. For a 
fraction of 21% of the total sample, the intrinsic [OIII] pseudo¬ 
equivalent width is in excess of 1 A, while about 1/7 of these 
galaxies (16 out of 498) display a striking (star-formation) activ¬ 
ity, with [OIII ]’5007 < — 2 A. 

ellipticals. Within the limits of lOgando’s et ahl (|2008h cal¬ 
ibration (see in particular Table 7 and 8 therein), in fact, 
these data indicate that a substantial fr acti o n of galaxies 
displays some residua l Oxygen emission llSarzi et alJ 12009 : 
IPapaderos et al.l 120131) . This important issue can be bet¬ 
ter assessed by means of Fig. |T] where we derive the dis¬ 
tribution of the intrinsic emission [OIII]' as from eq. (1311) 
The whole sample of 498 galaxies with available Fe5015 and 
[OIII ]5007 indices leads to a mean index strength {[OIII]') = 
—0.70 ± 0.54 A, with an evident skewness toward enhanced 
[OIII] emissions. Just as a reference, 16 objects out of 498 
(3%) in the figure display an intrinsic [OIII] feature stronger 
than 2 A, while for 107 galaxies ( 21 %) the [OIII] emission 
exceeds 1 A. 

After disentangling the combined behaviour of [OIII] 
emission and Fe5015 absorption, we are now in a better 
position to tune up an H /3 correction scheme. Again, our 
A-A diagnostic plot may add useful hints to tackle the 
problem. The rationale, in this case, is that the observed 
A [OIII] 5007 should directly trace the corresponding [OIII]’ 
intrinsic change, via eq. ED), providing the Fe5015’ strength 
is the same at symmetric distances from the galaxy center. 
We could therefore probe the A H/3 data set for any possible 
correlation with the corresponding A [OIII] 5007 output in a 
fully empirical way. 

For this task, we divided the whole galaxy set of Tabic [2] 
(Sample “A”) in different subsamples according to the “ac¬ 
tivity level” as traced by the overall spectral characteristics. 
In particular, we considered three refence samples. Sample 
“B” consists of the 7 “quiescent” galaxies with no reported 
emission lines. Sample “C” includes, on the contrary, the 
8 “moderately active” objects, marked with one asterisk in 
the table, wich are reportedly emitting at least the [Oil] 
3727 A feature. Finally, the 5 “most active” ellipticals (that 
is those marked with double asterisk in Table [2) are included 
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Table 3. Core [OIII]’ emission strength and observed AH(3 versus 
A [OilI] fitting slope 


Sample 

<[OIII]’)t 

a 

nt 

A - All galaxies (20 obj) 

-1.1 

0.08 

(68) 


±1.7 

±11 


B - Quiescent only (7 obj) 

-0.15 

-0.23 

(23) 


± 74 

±9 


C - Interm, active (8 obj) 

-0.67 

0.53 

(25) 


± 49 

±17 


D - Most active (5 obj) 

-3.1 

0.74 

(17) 


±2.3 

±46 



t Mean core [OIII]’ index, in A, as from Buzzoni et ah (2015) 
t Number of available pairs of points in the a fit 


in Sample “D”. They display [Oil] 3727 A emission accom¬ 
panied by the strongest [OIII]’ index in Fig. [S] Our results 
are summarized in the four panels of Fig. [8] As usual, the 
fitting slope after 3a data clipping is presented in Table [3] 
together with the mean of ce ntral corrected [OIII] emission, 
(as from iBuzzoni et alJl2015h . 

A firm important conclusion of our analysis, is that no 
striking correlation seems in place across the whole galaxy 
sample between H/3 and [OIII] 5 oo 7 - A mild positive slope is 
found overall (Sample “A”, upper left panel in Fig. [8]). only 
marginally in excess to the flat-slope case (a = 0.0-0.2). 
The flattening effect is clearly induced by the fraction of 
“quiescent” ellipticals (Sample “B”, upper right panel), for 
which we even notice a slight anti-correlation (a ~ —0.2). 
However, with increasing stellar activity, as in Sample “C” 
and “D” galaxies, one has to notice some hints of a change. 
Although with increased scatter and a much blurred trend, 
the data seem to point to a steeper slope a and, at least 
for Sample “D” (lower left panel in the figure), our results 

nominally match [Gonzalezj l) 19931 ) value. _ 

_ As a plain “thumb rule”, we could state that iGonzalePI 

dl993T) correction scheme fully holds for the most active 
strong-lined ellipticals, where [OIII]' emission roughly ex¬ 
ceeds 2 A. For intermediate cases (say for [OIII]’ emission 
strength about 1 A) a flatter slope (a ~ 0.3-0.5) might be 
more appropriate, while a negligible H/3 correction (a ~ 0.1) 
is required for weaker [OIII] emission. As shown in Fig. [9j 
a logistic curve could be adopted as a convenient analytical 
fit of this scheme, such as 


a = 


0.75 
1 + e k ’ 


(33) 


with k = ([OIII]' + l)/0.25, and [OIII] emission assumes, 
of course, negative values. The corrected H/3 index therefore 
derives as 


H/3coti = H/3 ohs - H/3 e = Hfiohs - a [OIII]'. (34) 


3.3 Mgi + Mgb = Mg 2 ? 

The prominent Magnesium feature about 5200 A is by far 
the most popular target in optical spectroscopy of elliptical 
galaxies an d oth er old stellar systems. As explored in full 
detail bv iMouldl (119781 4 . from a theoretical point of view, 
this feature is in fact a blend of both the atomic contribu¬ 
tion of the Mgb triplet at 5178 A, and the (1,1) A 2 n A' 2 E + 
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Figure 8. The A-A correlation test, as for Fi g. [3j for the 
[OIII] 5007 and H/3 indices in the adopted IBuzzoni et alj {2009b, 
Imi) sample of early-type galaxies. The whole Sample A, of the 
20 objects (upper panel), is split into three groups (i.e. Sam¬ 
ples “B”, “C” and “D”, from top to bottom) consisting, respec¬ 
tively of 7, 8, and 5 galaxies with increasing gas emission activity. 
Groups have been assembled according to the asteri sk ranking 
of Table [l] that combines the iBettoni &; Busoni dl987T) empirical 
classification with the [OIII] index streng th in our observa tions. 
The dashed line, in each plot, marks the lGonzale3 dl993h A-A 
correction slope. Our least-square fit is reported top right in each 
panel, and displayed (solid line) together with its corresponding 
Tier uncertainty (yellow fan). Our output is summarized in Ta¬ 
ble [3] 
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Figure 9. The devised correction scheme for H/3 in the form of 
a logistic curve versus intrinsic [OIII]’ emission, as from eq. (1331) . 
The “quiescent”, “active” and “most active” galaxy samples of 
Fig M are sketched by the three shaded boxes, with [OIII]’ mean 
central values and a error bars as from Table [3] 


vibrational band head of the MgH molecule. This transi¬ 
tion of Magnesium Hydrid e has a dissociation energy of 
1.34 eV (ICoelho et al.ll2005l . as from lHuber fe Herzberelll979l 
data), a much lower threshold than the first-ionization po¬ 
tential of 7.65 eV for atomic Mg. As a consequence, the 
Mgb absorption tends to be the prevailing one in warmer 
(type F-G) stars, while it accompanies the surging MgH 
absorption at the much cooler temperature regime of K- 
M stars. In addition, the relatively more flimsy structure 
of the molec ule m a kes MgH very sensitive to stellar lumi¬ 
nosity class (jBuzzoni et all [20011 1. as gravity acts on elec¬ 
tron pressure of the stellar plasma and slightly modulate 
the dissociation energy threshold. This dependence, how¬ 
ever, is severely tackled by the intervening appearence of 
TiO bands among M-type stars, which substantia l ly af- 


feet the pseudo-continuum level ( 

Tantalo & Chios! 20041: 

Buzzoni. Gariboldi, & Manteeazza 

19921, see, for instance 


Fig. 5 therein). 


The Magnesium complex is sampled by three Lick in¬ 
dices, namely Mgi, Mg 2 , and Mgb. The strategy is for 
Mg 2 to somewhat bridge the other two, with Mgi espe¬ 
cially suited to probe the molecular contribution, and the 
Mgb index better tuned on the atomic triplet (see Fig. 1101) . 
Within this configuration, Mg 2 shares its feature window 
with Mgb, and its pseudo-continuum windows with the Mgi 
index. From its side, Mgb aims at recovering the genuine 
atomic contribution by setting its pseudo-continuum level 
about the Mgi bottom (see again Fig. [TO] for an immediate 
sketch). 

Similarly to what discussed in previous section, we can 
further expand our scheme to assess in a more straightfor¬ 
ward way the Mg 2 -Mgb-Mgi index entanglement. First of 
all, as both Mg 2 and Mgi are expressed in magnitudes and 



Wavelength [A] 

Figure 10. An illustrative sketch of the Mg index set, and the 
two Fe5270 and Fe5335 indices. For each index, the blue and red 
sidebands are marked, flanking the relevant feature window. The 
nested configuration of the Mgb, Fe5335 and Mg 2 indices is well 
evident, as discussed in the text. The repr esentative spectrum of 
the elliptical galaxy NGC 4649 (from the iBuzzoni et ahll2009b[ 
data) is superposed for the sake of comparison. 


Mgb in A, it is convenient to recall the basic relations of 
eq. m and © and write 


' h/U 

— jQ-0.4Mg 2 

fl/fc 

— IQ-0.4 Mg x 

. fb/fl 

= (1 - Mgb/iUb) 


We can also combine the last two equations above to obtain 

fb = (fb\ fh\ Mgb 

A V/i ) \fc) V w b 

Quantities / 2 ,/i, and fb in previous equations refer to the 
flux density of Mg 2 , Mgi and Mgb feature windows, respec¬ 
tively, and Wb is the width of the Mgb feature window. For 
Mgi and Mg 2 this flux is normalized to the same pseudo¬ 
continuum, f c , while the fi bottom level can also be seen as 
a proxy of the Mgb pseudo-continuum. 

Likewise eq. m , also for the Mg case we can write 

/2 W 2 = fl {W 2 - w b ) + fi w b , (37) 

being W 2 = 42.5 A and Wb = 32.5 A (from Table © the 
width of the Mg 2 and Mgb feature windows, respectively. By 
matching eq. ([35). ([36]) and m and taking the logarithm, 
after little arithmetic, we lead to the final form that links 
the three indices: 

Mg 2 - M gl =-2.5 log (l - +k. (38) 

In the equation we left a fine-tuning offset ( k ) in order 
to account for the little aproximation we made in eq. ([35]) for 
fi to be an exact proxy of the Mgb pseudo-continuum level. 
It is interesting to verify, on empirical basis, the invariancy 
of the function 


^ l 0 -°- 4M si. (36) 
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Figure 11. An empirical check of the expected relationship 
among the three popular Mg indices sampled by the Lick sys¬ 
tem. The synthesis function <f)(Mg) of eq. (139t) has been scruti- 
nized by means of the two extended early-type galaxy samples of 
iTrager et al.1 J 19981 (406 galaxi es with complete Mg i, Mg 2 and 
Mgb data, see upper panel), and lOgando et al.l i 20081) (481 galax¬ 
ies, lower panel). The match with both samples allowed us to ob¬ 
tain the fine-tuning offset in the theoretical relationship, such as 
k = 0.03 mag. After offset application, notice for both plots the 
lack of any drift in the point distribution, a feature that assures 
the correct parameterization of our theoretical relationship. In 
addition, the point scatter around the zero line indicates that in¬ 
dex mutual relationship can be obtained within a ~1.5% internal 
accuracy. 


0(Mg) = Mg 2 - Mg, + 2.5 log 1 - 


Mgb 

W 2 


— k. 


(39) 


Again, the Ogando et al.l d2008l f together with the 
iTrager et al.l ([19981 ) galaxy database provide us with an ex¬ 
cellent reference tool. Our results are summarized in the 
two panels of Fig. [TT] As predicted, the lack of any drift 
of 0(Mg) versus Mg 2 (and versus Mgi and Mgb, as well) 
indicates that no residual dependency has been left unac¬ 
counted in our analysis. By searching the least-square fit 
to the data, we also get an empirical estimate for the fine- 
tuning offset, that is k = 0.03 mag. Providing to observe 
two indices, eq. (13811 allows us to secure the third one within 
a ~ 1.5% internal accuracy^ 



Figure 12. The A-A correlation test, as for Fig. [3] for the Mg in- 
dex complex. The AMqi, A A/c 2 and A Mgb variation across the 
iBuzzoni et al.l (l2009bl . 120 1 4) galaxy sample of Table [2] is probed 
in terms of the expected relationsip, as from eq. d40D (dashed line 
in the plot). Our least-square fit is reported top right, and dis¬ 
played (solid line) together with its corresponding ±lcr statistical 
uncertainty (yellow fan). 


A further check can be carried out in terms of our A- 
A test. By differentiating eq. ([38]). and considering that, in 
general, Mgb < W 2 , we obtain, in fact 

m S2 -«g,^ (1 ” 0.024 rfMgb, (40) 

The thoretical relationship is probed in Fig. [12] with our 
galaxy sample of Table [5] The usual least-square fit proce¬ 
dure to a total of 66 points provides a value of the slope 
a = 0.019±o.oo3- 


3.3.1 The /AT /5199 hidden bias 

As for the case of the [OIIIJsooz perturbing emission to 
the Fe5015 index, also for the Mgb index one may ex¬ 
pect a shifty, and mostly neglected, bias from overlapping 
gaseous emission lines. This is the case of the [NI] forbid¬ 
den doublet at 5198/5200 A (hereafter [Nljsigg), a some¬ 
what ominous weak feat ure that characte rizes the spec¬ 
tra of LINER galaxies dOsterbrockl 1 19741 . see Appendix 
4 therein), including e llipticals, either due to residual 
star formation activity (|Sarzi et al.l 120061 : ICapnellari et al.1 
(m3), or simply supplied by the post-AGB stellar compo¬ 
nent |Kalerl 198fll: Ferland et al.l 120121 : iKehrie et al.l 120121 : 
IPanader os et alJ l2013h. 

iGoudfrooii fe Emselleml d 19961 ) have been discussing in 
some detail the perturbing effect of [NIJsigg on the com¬ 
puted Mgb index in the spectra of elliptical galaxies. Differ¬ 
ently from the [OIII] 5 oo 7 versus Fe5015 interplay, the [NIJsigg 
emission positively correlates with the affected Mgb as [NI] 
luminosity enhances the pseudo-continuum level in the red 
side band of Mgb (see Fig. [TO]). thus leading to a stronger 
index. 

A simple correction scheme can be devised, to size up 


2 Note that, for vanishing values of Mgb as in low-metallicity 
stellar populations, eq. 1381) further simplifies as its r.h. term can 
be usefully aproximated by the linear term of its series expansion, 


i.e. [k-\- 1.09 (Mgb/u; 2 )] ~ 0.03 (Mgb +1). In this case, Mg indices 
straightforwardly relate as Mg 2 ~ Mg-L +0.03 (Mgb + 1). 
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Figure 13. The i nferred distribution of the “hidden” [NIJsigg 
gas emission in the lQgando et all d20Q8h galaxy sample. In case of 
intervening [NI] emission, the Mgb index is expected to enhance, 
on average, by the amount of eq. 03 with respect to its theoreti¬ 
cal figure, as from eq. J38l ). This excess (in the sense “Observed - 
Computed”, see lower x axis), and the related [NI] strength (up¬ 
per x axis) is probed by contrasting the subsample of 41 “active” 
galaxies (as traced by explicit [OIII] emission stronger than 1.5 A, 
as in the upper histogram) with the population of 275 “quiescent” 
galaxies with vanishing [OIII] (see the lower histogram). See text 
for a discussion. 

the [NIJsigg bias just relying on the third relation of eq. (1551) 
(and, again, assuming /i as the reference Mgb pseudo- 
continuum flux density). By passing to logarithms we have: 

hTi— = ln(/ 6 )-l n (/ 1 ) (41) 

V Wb J 

Again, considering that Mgb <C Wb, the l.h. term of the 
equation can be approximated by — (Mgb/ws) so that, by 
differentiating with respect to fi, we have 

rfM gb ^ <fh / 42 x 

w b fi ' 

If [NI] 5 i 99 emission appears, then a supplementary flux Fjv is 
added to the Mgb pseudo-continuum, and its reference flux 
density changes by d/i = F N /Wb, where Wb = 33.7 A, as 
from Table [7] As, by definition, the [NI] equivalent widtqjis 

3 Note that we maintain a negative notation for emission lines. 


[A 7 ] 5 i 99 ~ —Fn// i, eq. (l42l) can eventually be re-arranged 
in its final form 

dMgb = -[AT] 5 i99 k -0.96 [NI] 5 im . (43) 

Our result nicely compares with iGoudfrooii fc Emselleml 
dl996tl . who led to a more crude analytical estimate of the 
effect su ch as dMgb ~ —1.1 [A T] 5199 . 

The lQgando et all (i2008T ) extended galaxy database al¬ 
lows us a further independennt assessment of the [NIJsigg 
bias among elliptical galaxies. In this regard, we set up a 
simple test aimed at comparing the observed Mgb index for 
the entire galaxy dataset with the corresponding “predicted” 
output as from eq. (I3S[) . Galaxies with intervening [NI] emis¬ 
sion are expected to display (on average) a stronger-than- 
predicted Mgb index, and one could envisage to detect this 
signature in the index (O-C) (i.e. “observed - computed”) 
residual distribution. 

In order to get a cleaner piece of information, we 
only restrained our analysis to the galaxy sub-sample with 
stronger and poorer [OIII]soo 7 emission, assuming [O III] 5007 
to b e also a confident proxy for [NIJsigg activity dOsterbrockl 
11974t l. According to the [OIII ]’5007 distribution of Fig. [3 
we therefore selected a sample of 41 “active” galaxies with 
[OIII ]’5007 < -1.5 A emission, to be compared with a 
group of 275 objects with average or less-than-average (i.e. 
[OIII ]’5007 > -0.7 A) emission. For each galaxy with ob¬ 
served Mg 2 and Mgi indices we then computed the corre¬ 
sponding Mgb, via eq. (1551) (by setting k = 0.03 therein), 
and compared it with the observed entry in the catalog. The 
output distribution for the “active” and “quiescent” galaxy 
samples is displayed in Fig. 1131 According to eq. (1431) . the 
(O-C) residuals have been directly converted into a fiducial 
[NIJsigg emission strength, as in the upper x-axis scale of the 
figure. Within the uncertainty of our theoretical procedure, 
just a glance to Fig. [13] confirms the lack of any [NI] emis¬ 
sion among “quiescent” galaxies (([ATJsigg) = +0.15±o.O3 A, 
on average), while a weak but significant (> 3cr) activ¬ 
ity has to be reported, on the contrary for “active” sys¬ 
tems, with ([AT] 5199 ) = — 0.29±o.n A. Interestingly enough, 
the [NI]/[OIII] ~ 0.3/1.5 observed ratio for the “active” 
sub-sample matches in the right sense the corresponding 
empirica l figur e for_ the LINER galaxies, as suggested by 
IGoudfrooii fc Emselleml dl996l . see Table 1 therein). 


3.4 Spectral duality: the case of Fe5335 

Due to wavelength oversampling, along the Lick-index se¬ 
quence it may even happen that individual spectral features 
are accounted twice being also included in some pseudo¬ 
continuum windows of other indices. One relevant case, in 
this sense, is the Fe5335 index, which is fully comprised in¬ 
side the (red) pseudo-continuum window of both the Mg 2 
and Mgi indices. In addition, it also entirely shares its own 
red pseudo-continuum with the blue side-window of Fe5270 
(see, again, Fig. 1101) . Evidently, this makes both the Fe5270 
and the Mg indices to somewhat correlate with the Fe5335 


This explains the minus sign in the relevant equations throughout 
our discussion. 
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Figure 14. The A-A correlation test, as for Fig. [3] for the Fe5270 
and Fe5335 indices. The low correlation between index variation, 
as predicted by eq. J48l ) and d5 1 1) (dashed lines in the two plots, 
with their corresponding Tier Bayesian uncertainties, as from the 
blue fans) is consistently resumed by our observing test on the 
galaxy sample of Table [2] (solid lines with their corresponding 
±lcr statistical uncertainties, as from the yellow fans). See text 
for a discussion. 


strength. To explore this effect let us start first with the 
simpler case of Fe5335-Fe5270 entanglement. 

According to Table [T] the Fe5270 pseudo-continuum is 
sampled along a total of 47.5 A (i.e. red+blue side-windows, 
VF 52 ), whose 11.3 A are shared with the Fe5335 red side¬ 
band (W 53 d ). The total side-windows for Fe5335 (i.e. W 53 = 
Wlf + VF 53 lue ) amount to 21.2 A, and both features are 
sampled within similar wavelength windows (11152 and 1053 ) 
of 40 A each. By recalling eq. © and passing to logarithm, 
for Fe5270 we have 

In fl - ' Fe527 ° i = ln(/ 52 ) - ln(/| 2 ). (44) 

\ W52 ) 

As we are dealing with a weak absorption line (so that 
Fe 5270/«i52 <C 1) the l.h. term of the equation can be ap¬ 
proximated by — (7 ? e5270/w52). By differentiating, 

( dFe 5270 d/ 5 c 2 d/52 (45) 

' ^52 /s2 /52 

A similar notation holds for Fe5335, too. Any change (for 
whatever reason) in the /| 3 flux density, would therefore re¬ 


verberate into a change in the Fe5270 and Fe5335 strengths 
as well, so that 

( dFe 5270 N dfj 2 df^Wg_ ( dFe5335 \( Wjf \ 

V ™52 ) / 5 C 2 ~ H 3 ^52 ~ V W53 / V W 52 > ’ 1 ’ 

or 

dFe 5270 = ( A|) dFe 5335 = 0.238 dFe 5335. (47) 

The chance for Fe5335 to entangle Fe5270 depends on the 
probability for any random change in the Fe5335 index 
bands to affect the wavelength range in common with the 
Fe5270 pseudo-continuum. This can be be estimated of the 
order of Ws 3 A /(W $3 + W53) ~ 11.3/61.2 = 18% being, there¬ 
fore, 82% the probability for no entanglment at all between 
the two indices. By accounting for these figures, the average 
expected relationship between Fe5270 and Fe5335 eventually 
results in 


<£Fe5270 = (0.238 x 0.18) + (0 x 0.82) dFe 5335 
= 0.043±o.o9i dFe 5335. 


(48) 


This value has to be compared with the A-A test as dis¬ 
played in the upper panel of Fig. 1141 From our data (66 
points in total) we obtain a = 0.059±o.o68- 

Quite importantly, note that eq. m cannot be sim¬ 
ply reversed, once considering the impact of the overlapping 
window on the Fe5335 strength (because W 53 ^ VF 52 ). By 
reversing the indices in eq. (|46l). in fact we have 


/dFe5335\ 

d /53 

d /52 W$3 d 

/dFe5270\ 

(W£f\ 

V W 53 ) 

fs 3 

n 2 was 

V W 52 ) 

V IT 53 / 

or 





dFe 5335 = 


dFe5270 = 0.533 dFe 5270. 

( 


V 21.2/ 





), (49) 


(50) 


This interaction has a probability of the order of 
W$ 3 d /(W 52 + W 52 ) ~ 11.3/87.5 = 13%. so that, similar to 
eq. (148(1 . the eventual relationship, by reversing the indices 
in a A-A test becomes 


dFe 5335 = (0.533 x 0.13) + (0 x 0.87) dFe 5270 

= 0.07± o .is dFe 5270. 


(51) 


Again, this relationship can be probed by means of our A- 
A test, as shown in the lower panel of Fig. 1141 suggesting 
Q = 0.19±o,22 ■ 

As far as the Mg 2 - (or Mgi-) Fe5335 entanglement 
is concerned, we can proceed in a similar way, depend¬ 
ing whether a change occurs in the feature or in the side¬ 
band windows of the Fel index. According to Table [I] the 
Mg 2 (or Mgi) pseudo-continuum is sampled along a total of 
W 2 = 127.5 A. Of these, 40 A and 21.2 A are shared, re¬ 
spectively, with the feature ( 1053 ) and the side bands (W 53 ) 
of the Fe5335 index. 

By recalling that Mg 2 (and Mgi) are expressed in mag¬ 
nitude, the equivalent of previous eq. 03, in case of an 
intervening change in the Fe5335 pseudo-continuum flux, 
becomes 


dMg2 ~f 

f 2 fs 3 V W 53 > V W 2 > 

which leads to 

dMg, = ( - — - ) dFe 5335 = +0.0042^5335. 

2 \40 x 127.5/ 


(52) 

(53) 
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Figure 15. The A-A correlation test, as for Fig. [3] for the 
Fe5335 and Mg indices. The observed relationship for Mgi and 
Mg 2 is averaged, to gain a cleaner view of the effect. In spite 
of the average low correlation, as predicted by eq. d56l> (dashed 
line) one has to notice the important Bayesian uncertainty (blue 
fan in the plot, tracing the drier uncertainty), which makes the 
bias effect of Mg-Fe5335 entanglement possibly of non-negligible 
impact on individual galaxies. 


^ i i i i | i i i i | i i i i | i i i i 

teo: AMg =-0.004 ±0006 AFe5335 
obs: AMg = 0.002 ±0002 AFe5335 



J_i_i_i_i_I_i_i_i_i_I_i_i_i_i_I_i_i_i_L 


This outcome has a probability of the order of W^/iW^z + 
wss) - 21.2/61.2 = 35%. 

On the other hand, any change in the Fe5335 feature 
window will affect the Mg 2 (Mgi) strength by 


dMg 2 


<VS 

d /53 w 53 

/dFe 5335 \ 

(W$ 3 \ 

/I 

fss W2 

V W53 ) 

\w 2 ) 


(54) 


By replacing the relevant quantities, in this case we have 
dMg 2 = - ( ) dFe 5335 = -0.0078 dFe 5335, (55) 


with a probability of ws3/(Ws3 + W 53 ) ~ 40/61.2 = 65%. 
Therefore, if a random process (of whatever origin) is at 
work by changing Fe5335, then the weighted average of 
eq. (1531) and (1551) eventually provides us with the expected 
relationship for Mg 2 (and similarly, for Mgi, as well): 


dMg 2 = dMgi = -0.0036±o.oo57 dFe 5335. (56) 

Our theoretical output is probed in Fig. 1151 by averag¬ 
ing the AMgi and AMg 2 versus AFe5335 variations for a 
cleaner display of the effect. Observations actually confirm 
the nearly vanishing relationship in place across the whole 
galaxy sample, making the expected impact of Fe5335 ran¬ 
dom change on Mg indices negligible. 

On the other hand, if a change in the Fe5335 feature 
occurs as a result of a genuine (selective) variation of Fe 
abundance, then only eq. (ESI) has to be taken as a reference. 
Although small, the induced bias on Mg 2 (and Mgi) may 
be an issue in this case, expecially if we aim at carefully 
investigating the [a/Fe] ratio in a galaxy sample. In fact, any 
stronger Fe5335 index would induce, by itself, a shallower 
Mg index. 


3.5 Mr. Iron, I suppose... 

About half of the total Lick indices are meant to trace Iron 
as a main elemental contributor to the index feature. The 


ominous presence of weak (i.e. non-saturated) absorption 
lines of Fel and Fell in the optical/UV wavelength region is, 
of course, a well established property of stellar spectra (e.g. 
|Pagel|[l997l h that can be exploited by high-res spectroscopy 
to set up the linear branch of the spectral curve of growth, 
thus leading to an accurate measure of the elemental abun¬ 
dance. This figure is usually taken as an empirical proxy 
of stellar “metallicity”, at large, assuming Iron to scale ac¬ 
cording to all the other metals. On the same line, also UBV 
broad-band photometry may rely on the “blanketing” effect, 
mainly driven b y the Iron con t ribution, to constrain stellar 
metallicity (e.g. lSandagelll969l : lGolavlll974l '). 

As far as low-res spectroscopy or narrow-band pho¬ 
tometry is concerned, as for the Lick indices, however, 
such a huge number of low-strength Fe lines may reveal a 
quite tricky feature, as Iron markers actually come heav¬ 
ily blended with other elements. The potential drawback 
has already been emphasized, at least in a qualitative way, 
in earlier works dea l ing with the Lick-system assessment 
dBurstein et al.lll984l : iFaber et al.lll985l : iGorgas et all 1 19931 : 
IWorthev et al. 119941 . amon g others). This analysi s has then 
been f u rther carried on bvlTripicco fc B ell (1995); Worthe^ 
(1996); iTrager et al.l (Il998tl , and Serven, Worthev. &; Brilev 
12005), urging in some cases a revision of the index label, as 


for the original Fe4668 index dWorthevlll996l : ITrager et al.l 
imi), now redefined as C 2 4668 according to the real pre¬ 
vailing contribution in the evolutionary context of old stellar 
populati ons, as in elliptical galaxies. 

The ITrager et al.l dl99af revision process has been in¬ 
tentionally left at the highest conservative level, likely 
to avoid any massive re-definition of the original index 
set. However, w e want to retake here some arguments of 
IWorthevAl dl996T ) preliminary analysis to show that, at least 
for some relevant cases dealing with bona fide Fe features, an 
important update effort seems yet to impose. We will defer 
this delicate task to a forthcoming paper discussing here, 
just as an illustrative excercise, the case of index Fe4531. 


3.5.1 Absolute and net index sensitivity 

The main rationale for Lick-index set up dWorthev et al.1 
mu) aims at catching all the prominent absorption fea¬ 
tures in the optical spectra of late-type (type G-K) stars. 
This uniquely requires the feature to be easily recognized 
at a ~ 8.5 A (FWHM) resolving power. If a blend could be 
envisaged for a feature, the latter is tacitly assumed to be 
named after the element contributing with the deepest line, 
as re cognized from the l ist of theoretica l atomic transitions 
(e.g. ISmith et al] 1 19961 : iPiskunov et al.l Il995l : iKupka et al.l 
Il999ll . The deepest line also defines the feat ure centroid, 
which is reported in the IWorthev et akl (Il994l l original Ta¬ 
ble 1. No special requirement is demanded to the two index 
side-bands, set to sample a nearby pseudo-continuum, other 
than being suitably “feature-free” and sufficiently extended 
in wavelength such as to be insensitive to stellar velocity 
dispersion broadening dWorthev et al.lll994lh 

So, while the main focus is on the feature prominence, 
that is on feature’s absolute sensitivity to the “leading” 
chemical element, one should not neglect the fact that it 
is the index “net” sensitivity (that is by netting the feature 
sensitivity to a given element with the impact of the latter 
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Figure 16. The Fe4531 index sensitivity to the 40 most prominent chemical elements (as sequenced in the x axis, disregarding 
their ionization states) in the Sun (left panel) and Arcturus (right panel) spectrum. In both panels, the yellow histograms refer to the 
elemental contribution (<S^) to the feature window, according to the definition of eq. JED, while the blue histograms are for the elemental 
contribution (<S^) to the flanking pseudo-continuum index sidebands, according to eq. J58D . The net contribution to the index strength 
(A<Sx), as from eq. <l59t) , is displayed in the lower plots, where “proactive” elements are displayed in red and “reactive” elements in grey. 
Note the leading role of Titanium and the reactive action of Vanadium and Magnesium. See text for a full discussion. 


on the pseudo-continuum bands, too) that eventually set the 
real index response to a given chemical species. 

Like lServen, Worthev, fc Brilev! (|2QQ5h , in order to ex¬ 
plore this important point we performed a test by prob¬ 
ing the spectral wavelength region around the Fe4531 fea¬ 
ture in the high-res spectra of two reference templates for 
old stellar populations, namely the Sun (type G2V, seen at 
R ~ 350, 000 resolving power) and the star Arcturus (type 
K1.5III, observed at R ~ 80,000). For our analysis we re¬ 
lied on the Spectro web 1.0 interactive data base of spectral 
standard star atlases llLobel2007l . l2008l . 1201 il l. The Web pro¬ 
gram interface allows an accurate fit of the observed spectra 
of these stars, supplying a detailed identification of all the 
chemical species with absorption features within the wave¬ 
length range (feature plus pseudo-continuum windows) of 
the Lick Fe4531 index. 

The sensitivity of both the feature- and the contin¬ 
uum windows to the first 40 most important absorbers has 
been quantified by summing up all the relevant line depths 
(d = 1 — 7(A)// C (A)) as from the high-res spectrum. As we 
always are dealing with weak (i.e. non-saturated) lines, one 
could assume, with sufficient accuracy, the equivalent width 
to simply scale with the line depth. So, a measure of the 
absorbing strength of element X (disregarding its different 
ionization states, in case) in the feature window of the index 
(of width w, as from Table [U can be written as 

W 

S f x = (57) 

The same can be done for the pseudo-continuum sensitivity 
to the corresponding element: 



where the width W comes again from Tablcjl] As, in general, 
w ^ W, in the latter case Sx should be normalized by a 
factor (w/W) to consistently compare with S* x , as displayed 
in eq. (ISHll . 

Clearly, the net sensitivity of the index to element X is 
probed by the difference 

A Sx=S f x ~S c x , 

or 

ASx d x J2 W d. 

w w W 

If ASx > 0, then element A' is “proactive” to the index 
strength (that is the index strengthens if A' abundance in¬ 
creases) otherwise, for A<Sx < 0, X is “reactive” (that is the 
index weakens if X abundance increases, as A more strongly 
absorbs the pseudo-continuum than the feature window). 
The results of our test on the Fe4531 index, are displayed in 
the two panels of Fig. 1161 for the Sun and Arcturus, respec¬ 
tively. 

Accordingly, in Table [4] we list the five most prominent 
absorbers to the feature window alone in the Sun and Arc¬ 
turus. These data definitely show that Fe is the strongest 
contributor to the 4531 A blend itself. On the other hand, 
Fig. H6] also show that Fe lines affect with a similar strength 
also the pseudo-continuum flux density. When netting both 
contributions, as in Table [5] one is left with a quite dif¬ 
ferent list of prevailing chemical species that really settle 
the Fe5431 behaviour. In particular, the main driver of the 
index appears to he Titanium (!), followed to a lesser ex¬ 
tent by Chromium. Note, in addition, the weak but steady 


(59) 

(60) 
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Table 4. The five chemical species most important contributors 
to the 4531 A blend in the Sun and Arcturus spectra 


Sun 


Chemical 

species 

St 

Chemical 

species 

st 

Fe 

22.3 

Fe 

44.0 

Ti 

13.5 

CN 

35.3 

Cr 

10.1 

Ti 

25.6 

Ni 

2.6 

Cr 

23.7 

Co 

2.2 

Ce 

11.2 


t According to eq. eg. 


inverse correlation of the index with the Vanadium abun¬ 
dance, which acts therefore as a “reactive” element. 

We could further carry on our analysis by probing the 
Fe4531 index response to homogeneous groups of chemical 
elements, like for instance the q-elem ent cha in, more di¬ 
rectly related to SNII nucleosynthesis llWooslev fc Weaver! 
1 19951 ). (including the combined contribution of Ca, Mg, 
MgH, Na, Si, SiH), or the Fe-Ni group (including Co, Cr, 
Fe, Mn, Ni, V), as a main by-product of SNIa activity 
dNomoto et al.l Il997f) . This is summarized in the bottom 
lines of Table [5] While Ti is definitely the main element 
that modulates Fe4531 strength, the other SNII a elements 
act as very weak “reactive” contributors, opposite to SNIa 
elemental group, which “proactively” contributes to the in¬ 
dex strength, mainly in force of Cr, Co and Ni contributions. 
Definitely, in spite of its outstanding role, our conclusion is 
that Fe plays as a quite dull actor in this game. 


4 DISCUSSION fc CONCLUSIONS 


Since its definitive settlement, in t he early 90ies 
ll Coreas et al.l fl993l : IWorthev et al.l IT 9941) . the Lick-index 
system has broadly imposed as the standard in the study of 
stars and old stellar systems. As a winning strategy, it linked 
the low-resolution spectroscopic approach to the narrow- 
band photometry such as to deliver at a time a condensed yet 
essential piece of information about all the relevant absorp¬ 
tion features that mark the SED of an astronomical target 
in the 4000-6000 A wavelength range. 

Nowadays, the Lick-index analysis consistently 
flanks other resources of modern extragalactic astron¬ 
omy at optical wavele ngth, including high-resolution 
spect r al synthesis (e.g. IVazdekisI |200l|^Buzzonie^_aL| 
2005^Rodriguez^Merino^^ahLJ^OOSL^ertonee^alJ 20081 : 


Vazdekis et al.l 20ld : iMaraston fc Strombackl 20111 1 and 


tribution (SED) (e.g- Massarotti, Iovino, & Buzzonil 

2001 ; 

Bolzonella et alJ 201$"" Dominguez Sanchez et al. 

2011 ; 

Pforr. Maraston. fc Ton ini 2012; Gonzalez et al. 

2012 ; 

Scoville et al. 20131). However, as we have being discussing 


before, none of these is free from some limitations so that 
narrow-band spectrophotometry still remains, in most 
cases, the only viable diagnostic tool to tackle with some 
confidence the evolutionary status of unresolved distant 
galaxies. 

Along the original bulk of 21 indices, later on 
(IWorthev fc Ottavianilfl997f) , the Lick system has been com¬ 


Table 5. The five chemical species with the most prominent net 
sensitivity in the the Fe4531 index for the Sun and Arcturus spec¬ 
tra 


Sun 

Arcturus 

Chemical 

species 

AS+ 

Chemical 

species 

ASt 

Ti 

7.9 

Ti 

9.8 

Cr 

3.3 

Cr 

4.0 

V 

-1.4 

V 

-2.7 

Ni 

1.1 

Co 

1.9 

Mg 

-0.9 

Sm 

1.8 


SNIa group 

4.4 

2.7 

SNII group 

-1.4 

-1.0 


I According to eq. 1591 ). 


plemented by including some higher-order Balmer lines 
(namely, Hy and H<5), with the aim of overcoming the possi¬ 
ble perturbation of gas emission on the other optical fea¬ 
tures (especially H/f) in the galaxy spectra. An alterna- 
tiv e approach to the same problem was already attempted 
by iGonzalej dl993l h to size up the H/3 emission compo¬ 
nent in the spectra of active elliptical galaxies, by using 
[OIII] 5007 A emission as a proxy. Although nominally 
not part of the original Lick system (being the only index 
that explicitely covers an emission feature), the lGonzale3 
d 19931) [OIII] index is often included in the extragalac- 


tic studies accompanying more classical Lick analysis 

(e.g. 

Kuntschner et al. 200ll; Prugniel. Maubon. & Simien 

2001 |; 

Beuing et al. 20(53; Nelan et al. 20051; Ogando et al. 

200 $ 

Wegner & Grogin 20081). To all extent, we can consider it as 


a further addition to the standard system, which amounts, 
allover, to a total of 26 indices (see their wavelength distri¬ 
bution in Fig. [T}. 

A recognized difficulty, when using Lick indices from 
different observing environments, deals with the standard¬ 
ization procedure. As we have being discussing in Sec. 2.1, 
such a trouble comes out, from one hand, from the inherent 
difficulty with the Lick IDS set of primary stellar calibrators, 
originally ge nerated by the analogical output of a video cam¬ 
era tube dRobinson fc Wamplen 19721 ). As a result, IDS spec¬ 
tra lack any stable wavelength scale and intensity response, 
a limitation that, by itself, constrains calibration accuracy 
and p revents any consistent assessment of the spe ctra S/N 
ratio dFaber fc JacksorJll976l : IWorthev et al.lll994 ). 

On the other hand, the problem of index standardiza¬ 
tion intimately relates to a tricky and yet not fully realized 
theoretical drawback. In fact, even if both feature (/) and 
pseudo-continuum (f c ) density fluxes, as of eq. ([l]), obey a 
normal distribution, the 7Z = (// fc) ratio may notQ The 
index definition, in terms of ratio TZ, is a direct legacy 
of the “spectroscopic mood” of the Lick system, and it 
closely recalls the standard procedure to derive line equiva¬ 
lent width as in high-resolution spectra. Given the different 
conditions of low-resolution unfluxed spectra, however, we 
have demonstrated that a somewhat alternate “photometric 
mood” should be preferred, with feature strength assessed in 


4 In some cases, the TZ r atio may even display a bi-modal statis¬ 
tical distribution, 1 s e e I M a r s i qf i allT 96 4| . Il965l) . 
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terms of ratio IZ' = (/ // c ). The latter is a far more suitable 
and robust statistical estimator, which displays a Gaussian 
distribution and turns index strength in terms of a narrow- 
band “color”. Under certain conditions, we have shown that 
7 Z tends to 7 Z 1 providing a spectrum is nearly “flat” and 
overcomes a minimum S/N threshold, such as 

/ a \ 1/2 

(S/N) oba > 22 (^—J ~50 1/2 \px~\ (61) 


being 0 the spectral dispersion (that is the observed wave¬ 
length pixel scale), and W the average index window size as 
from eq. (1211) and Tablc[l] This means, for example, that for 
a resolving power R ~ 2000, a minimum ( S/N) obs ~ 5 px _1 
ratio has to be reached, assuming a Bayesian pixel fair sam¬ 
pling. 


A relevant conclusion of our analysis is therefore that, 
in an effort to revise the Lick system, one should better 
consider to compute indices either in the classical form of 
pseudo-equivalent width, relying on the classical IZ ratio of 
eq. m, but from previously linearized spectra according to 
eq. m, or (more safely) in the form of “narrow-band” col¬ 
ors of simply fluxed spectra, that is in terms of the IZ' ratio 
of eq. m, and express them in magnitude scale, as from 
eq. J 2 |). 


As a further issue in our analysis, in Sec. 3 we assessed 
the possible redundancy in the index delivered information 
when matching for instance galaxy spectra. About 25% of 
the covered wavelength range, in fact, contributes to the 
definition of two or more Lick indices, a feature that may 
lead index changes to correlate, in some cases, beyond any 
strictly physical relationship. A relevant case has been ex¬ 
plored, for instance, dealing with the nested configuration 
of the Fe5335, Fe5270 and Mg 2 indices, a classical triad of¬ 
ten involved in literature studies to define “global” metal¬ 
lic ity ind ices, su c h as th e (Fe) = [Fe _5270 + Fe 5335)/2 
of lGorgas, Efstathiou, fe Aragon Salamancal (Il990h . or the 
perused \MqFe 1 = [ Mgb (Fe5270 + Fe5335)/2 ] 1 / 2 meta¬ 
index of Gonzalez il993j) , in i ts man ifold variants (see, e.g. 
Fritze-v. Alvenslebenl 1 19981: iThomas. Maraston. fe Bender! 


20031 : Clemens et al.l 2006ll . We have shown, in Sec. 3.4, that 


any intervening change in the Fe5335 strength also reflects 
in the Mg 2 index (and in Mgi, as well) in quite a puzzling 
way. From one hand, if a flux addition increases the pseudo¬ 
continuum level in common, then both Mg 2 and Fe5335 
strengths are seen to increase, according to eq. (1551) . On the 
other hand, if the Fe5335 index becomes stronger due to a 
deeper corresponding feature, then Mg 2 is seen to decrease, 
as from eq. dig). 

To further embroil the situation, in Sec. 3.3 we made a 
point dealing with the the Mgb index as a tracer of the 
galaxy metallicity, fo r instance wi thin the (MgFel meta¬ 
index. As shown by iGoudfrooii fc Emselleml (ll996l l. the 
strength of the Mg atomic feature itself could be tacitely 
modulated by the intervening gas emission, even in rela¬ 
tively quiescent stellar environments. The forbidden [NI] line 
emission about 5199 A adds flux to the red Mgb pseudo¬ 
continuum (see the sketch in Fig. 1101) . thus leading to a 
stronger Mg absorption feature, overall. In Fig.[l3]the effect 
has be en c learly detected on a statistical basis by relying 
on the lQgando et aD (120081 ') extended galaxy database. The 
[NI] emission is seen to correlate with the [OIII] strength, 
roughly in a ratio 1:5, and it appears to be in excess of about 


0.4 A among the most “active” (i.e. [OIII]-strong) objects in 
the sample compared with the quiescent ellipticals. Accord¬ 
ing to eq. 63. the Mgb strength is inc reased by a si milar 
figure . When relying, for instance, on the lPipino fe Danzigeil 
((2011) empirical calibration to convert Lick indices into 
metal abundance (see, in particular Table 5, therein), this 
leads to predict a slightly higher metallicity, even in excess 
of some 0.2 dex in star-forming galaxy environments. 

Like the [NI] case, hidden emission can subtly plague 
the nominal diagnostic scheme for other Lick indices. In¬ 
deed, the most debated issue, in this regard, deals with a 
fair correction of H/3 absorption in the spectra of ellipti¬ 
cal galaxies. Compared to all the other Lick indices, H/3 
(and its weaker B alrn er co mpanions, like H 7 and H<5, see 
IWorthev fe Ottavianil Il997f ~) is the only feature that pos¬ 
itively correlates with the effective tem perature of stars 
(see, e.g. Fig. 4 in iTrager et al.l Il998l h As it becomes 
strongest among A-type stars, this makes the index of 
special interest to break model degeneracy in the age di- 
agn ostic of old- and intermediate- (t <1 Gyr) age galax- 


ies (Buzzoni, Mante^azza, & Gariboldi 19941; 

Worthcy 1994); 

Maraston & Thomas 20001; Tantalo & Cliiosi 

2004). On the 

other hand, as gas and stars within a galaxy 

are conflicting 


in their line contribution (the first, by emitting, the latter 
ones by absorbing), the integrated H/3 strength displays the 
net effect of this interplay. _ _ 

In Sec. 3.2 we tackled the lGonzale3 (Il993li well expe¬ 
rienced correction procedure, that relies on the measure of 
the adjacent [OIII] forbidden emission at 5007 A, as a proxy 
to infer the H/3 gaseous contribution. Several physical argu¬ 
ments can be advocated to question such a straightforward 
relationship between the two emission lines, in consequence 
of a range of different evolutionary mechanisms (in primis 
LINER activity) at work. However, one can admittedly as¬ 
sume [OIII] emission to be a generic (but not exclusive) 
tracer of star-forming activity (and the implied presence of 
HII gas clouds) within a young stellar population. 

Taking a set of spectroscopic observations of 20 
bright elliptical galaxies, originally aime d at stu dyin g ra¬ 
dial spectrophotome t ric g radients (ICarrasco et al.1 1 19951 : 
iBuzzoni et al.l l2009bl . 120151 ), we devised a simple and ro¬ 
bust diagnostic scheme to probe the possible [OIII] versus 
H/3 entanglement. As observations were taken symmetri¬ 
cally across the galaxy major axis, one can reasonably as¬ 
sume that, whatever the distinctive evolutionary properties, 
the same stellar population (and, on average, gas content) 
is sampled at the same radial distance, on opposite sides 
from the galaxy center. If this is the case, then the H/3 and 
[OIII] 5007 index differences between the two galaxy sides can 
be contrasted in search for possible statistically correlations. 
One major advantage of this approach is that it is nominally 
free from zero-point and other calibration problems with the 
data, just relying on the assumed central symmetry of the 
galaxy stellar population. 

Our results, summariz ed in Fig. 0 and Table [3 
basically confirm I Gonzalez\ (l99&) correction figure of 
AH/3/A[OIII]5 oo7 ~ 0.7 ± 0.5, but only for those ellipti¬ 
cals with established on-going star-formation activity and 
stronger core emission (i.e. [OIII] 5007 /$ — 1.5 A). As far 
as more quiescent galaxies are considered, any H/3 emission 
is found to fade out more quickly than [OIII] 5007 emission, 
and this suggests a milder correction scheme, that we pa- 
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Figure 17. The devised scheme for H ft correction, as in Fig. [9j is applied to the lOgando et alj d2Q08T) early-type galaxy sample, 
displayed in the H/3 versus Mg 2 index plane. The strength of H (3 e gas emission is sized up and subtracted according to eq. 03. taking 
the observed [OIII ]’5007 emission as a proxy. To better appreciate the correction effect, we adopted three different galaxy markers for 
nearly “quiescent” (i.e. [OIIIJ ’5007 > — 1 A, or H/3 e > — 0.4 A, solid dots), “moderately active” (—2 <[OIII ]’5007 < — 1 A, or —1.5< 
H (3 e < — 0.4 A, triangles), and “active” (i.e. [OIH l ’soor < —2 A or H/3 e < — 1.5 A, rombs) galaxies. As a guideline, in both panels we also 
displayed the lBuzzoni, Mantegazza, &; Gariboldil i 19941 ) synthesis models for a standard Salpeter IMF and [Fe/H] metallicity parameter 
between —0.25 (upper line) and +0.22 dex (lower line). “Ladders” across model sequences tick galaxy age for values of 2, 5, 10, and 15 
Gyr, with increasing Mg 2 - After correction (right panel), galaxy H/3 index is increased, on average, by 0.25 d= 0.20 A along the whole 
sample, making galaxies to appear, on average, about 20-30% younger and 0.2-0.3 dex metal poorer. 


rameterized by means of a logistic curve (see Fig. [9j, as in 
eq. ([33]) and m- 

The effect of our H 0 correction c an be appreciated 
by looking again at the lOgando et alj J2008h sample of 
elliptical galaxies, as in Fig. 1171 The corrected galaxy 
distribution in the H/3 versus Mg 2 plane is shown in 
the right panel of the figure, compared with the origi- 
nal sample, as in t he left panel. J ust a s a guideline, the 
iBuzzoni. Mantegazza. fe Gariboldl d 19941 1 models are over¬ 
plotted, throughout. One can notice that most of those 
galaxies with “shallow” (or negative) H j3 in the left panel 
are effectively recovered by the procedure, leading eventu¬ 
ally to a much cleaner H/3 versus Mg 2 relationship. Along 
the whole galaxy sample, the corrected H/3 index slightly 
increases, on average, by AH/3 ~ 0.25 ± 0.20 A. Although 
apparently negligible, this figure has important impacts on 
the inferred galaxy age, which typically shifts toward 20-30% 
younger values with a 0.2-0. 3 dex poorer inferred metallic ity 
I Buzzoni. Mantegazza. fe Gariboldilll994l : IWorthevlll994l l. 

Aside the H/3 problem, the striking importance of the 
[OIII] 5 oo 7 forbidden line also emerged along our discussion 
on Sec. 3.2 when dealing with the intrinsic strength of the 
Fe5015 index. From one hand, any intervening [OIII] 5 oo 7 
emission will decrease Fe5015 strength as it adds flux to 
the feature window of the latter. On the other hand, the 
broader and stronger Fe5015 absorption in the integrated 
spectrum of old stellar populations may artificially depress 
the apparent [OIII] 5 oo 7 luminosity, thus leading to underes¬ 
timate its equivalent width. In our analysis, we tackled in 
some detail the problem of a mutual correction of the two 
indices devising a general procedure that could straightfor¬ 
wardly be applied to other similar contexts, as well. The in¬ 
trinsic [OIII] 5 oo 7 strength can therefore be derived in terms 


of observed quantities, via eq. m, while the intrinsic Fe5015 
index, corrected for O-emission, simply derives from eq. m- 

One important consequence of our analysis is that cor¬ 
rected Fe5015 index is generally expected to increase in 
the galaxy spectrum, opposite to Mgb trend. This feature 
may have some impact when combining the two indices, 
like in the [MgFe50] meta-index, so extensively used for 
instance in the Sauron studies dKuntschner et all 120101 : 
iFalcon-Barroso et al.l 1201 il l . A n intersting compari son can 
be done, in this regard, with [Peletier et ali (|2007l . see, in 
particular, Fig. 3, therein). In the plot, when taking the 
iThomas. Maraston. fc Bender] (l2003lf stellar population syn¬ 
thesis models as a reference in the Fe5015 vs. Mgb do¬ 
main, it appears that Sauron ellipticals better match the 
a-enhanced models with respect to the standard solar chem¬ 
ical partition. On the other hand, one may argue that this 
is the natural consequence of a relatively “strong” Mgb and 
“shallow” Fe5015 so that, when accounting for the envis¬ 
aged corrections, the net effect on the fit is to mitigate any 
inferred a-element enhancement of the galaxy population. 

For its relevance, the Mg index complex (including 
Mgb, Mgi, and Mg 2 indices) has been the focus of our 
specific analysis in Sec. 3.3. In its original intention, the 
Lick system aimed at singling out the atomic and molec- 
ular contri b ution to the bro ad “M g valley” about 5150 A 
dFaber, Burstein. fc Dressier] 1 1977h . This evidently leads to 
conceive some entanglement among the three indices, that 
we theoretically assessed, for the first time, by means of 
eq. (13811 . This inherent relationship can be usefully relied 
upon in case of incomplete spectral observations, either to 
derive the missing index (within a ~1.5% internal accuracy) 
given the other two or, even more importantly, to assess any 
possible deviation from the standard scheme, as due to in- 
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tervening bias effects in the spectra. This has actually been 
our approach, for instance, to single out the [NIJsigg contri¬ 
bution, as shown in Fig. 1131 

As far as low-resolution spectra are dealt with, as in 
any Lick analysis, we are left with the unescapable conse¬ 
quence that any absorption feature is in fact a blend of a 
coarse elemental contribution. As we have shown in Sec. 
3.5, this problem has been widely explored, from a theo¬ 
retical point of view, in th e rec ent literature (the work of 
IServen, Worthev. fe Brilevll2005l . is probably the most up-to- 
date attempt in this sense) in order to assess index “respon¬ 
siveness” to the different chemical species that may poten¬ 
tially intervene at the relevant wavelength of the blended ab¬ 
sorption feature. On the other hand, one has to notice in this 
regard that, for its inherent definition, Lick-index strength 
directly stems from the residual “in-band” elemental con¬ 
tribution to the feature itself, compensated by the corre¬ 
sponding “off-band” contribution to the flanking pseudo¬ 
continuum. Contrary to high-resolution spectroscopy, there¬ 
fore, a tricky situation may actually set in, where the Lick- 
index christening element (that is the one with the deepest 
absorption within the feature window) is not necessarily the 
representative one that constrains the index strength, as a 

whole. _ _ 

As a recognized e xamp le dTrager et alJ 1 19981 : 
IServen, Worthev. fe Brilevl l2005li . in Sec. 3.5 we re¬ 
took in somewhat finer detail the case of Fe4531, one among 
the many bona fide Iron tracers in the optical spectrum of 
galaxies, according to the classical Lick scheme. For our 
analysis we relied on the Sun and Arcturus observed spectra 
as typical main-sequence and red-giant stellar templates, in 
an attempt to extend our conclusions also to the old-galaxy 
framework, according to more elaborated stellar population 
synthesis models. The comparison of Table [5] and [5] is 
illuminating to catch the sense of our analysis. 

One sees, in fact (Tabic [4j that, both in the Sun and 
Arcturus (and, by inference, in the spectrum of elliptical 
galaxies, for instance), Iron is undoubtedly the leading con¬ 
tributor to the absorption feature at 4531 A, but once clear¬ 
ing the Fe coarse absorption in the nearby pseudo-continuum 
as well, it is actually Titanium that, almost solely, com¬ 
mands the game (Tabic [5j and eventually modulates Fe4531 
index behaviour. As a final consideration, the case of Fe4531 
and more generally the essence of all our previous discus¬ 
sion, evidently urge an important rethinking process of the 
Lick-index tool, to bring it fresh life and even better tuned 
diagnostic performance in view of its updated application to 
the new-generation spectral data to come. 
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